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SEISMIC  COMMUNICATION  IN  BASIN  AND  RANGE  PROVINCE  VALLEYS 


1.  INTRODUCTION 

1 . 1  Statement  of  Need 

"“It  has  long  been  recognized  that  the  Earth  provides  a  virtually  Inde¬ 
structible  full  duplex  path  for  communicating  with  hard  burled  facilities  <1). 
Work  on  the  subject  has  highlighted  estimates  of  channel  capacity  for 
communicating  In  competent  rock  (2),  as  well  as  the  performance  of  specific 
hardware  and  coding  at  short  and  long  ranges  (3»4). 

There  Is  a  need  to  determine  the  Information  capacity  for  communicating 
selsmlcally  In  a  landform  suitable  for  basing  MX  (5).  Of  paramount  Interest 
Is  the  Information  rate  for  communicating  selsmlcally  In  deep,  dry  alluvial 
valleys  of  the  Basin  and  Range  Province,  the  leading  landform  for  MX  basing. 

To  this  end  we  establish  the  maximum  Information  rate  for  sites  In  Railroad 
Valley,  Steptoe  Valley  and  Jackass  Flats,  Nevada  over  a  range  of  source 
strengths,  bandwldths . and  distances.^ 

1.2.  Approach 

We  treat  the  Earth  to  be  a  linear,  time  Invariant  transmission  path 
defined  by  Its  Green's  Function.  Seismic  signals  observed  at  a  distant 
receiver  point  are  delayed,  colored  versions  of  those  generated  by  the  source 
masked  partly  by  additive  noise.  For  surface  sources  and  surface  observers, 
seismic  transmission  at  a  distance  Is  dominated  by  waves  attenuated  by  Q 
losses  and  cylindrical  spreading. 

Seismic  noise  at  the  candidate  sites  Is  well  represented  as  stationary. 
Independent,  Gaussian  processes  convolved  with  the  appropriate  site  sensitive 
shaping  operators. 

Maximum  capacity  seismic  channels  are  computed  after  Hartley-Shann'''\ 
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for  the  seismic  signals  and  noise  terms  found  in  each  area.  At  high  signal- 
to-nolse  ratios  the  optimum  transmission  channel  is  tantamount  to  matched 
filtering.  The  effect  of  source  and  receiver  arrays  on  communication  is  also 
considered  for  the  simple  case  of  coherent  signals  and  incoherent  additive 
noise. 

2.  FINDINGS 

Coitmunlcatlon  by  seismic  means  is  as  commonplace  as  a  “knock  at  the 
door";  its  feasibility  is  not  in  question.  What  is  in  question  is  the  capaci¬ 
ty  of  seismic  data  links  operated  at  specific  source  strengths  and  distances 
in  the  Basin  and  Range  Province  Valleys. 

Information  capacity  for  a  ground  path  is  readily  calculated  once  the 
signal  and  noise  spectra  are  specified.  Data  links  can  be  optimized  by 
weighting  at  the  source  and/or  receiver.  In  any  event,  signal  distortion  due 
to  seismic  propagation  can  be  nullified  by  a  linear  operator  that  is  the 
inverse  of  the  seismic  impulse  response,  G(r,t). 

Seismic  waves  excited  by  surface  impacts  are  well  represented  by 
unimodal  Q  damped  surface  waves.  Under  this  representation  signal  spectra  can 
be  generated  over  any  des:.*  u  range  of  distances,  source  strengths  and 
bandwldths. 

The  effect  of  range,  bandwidth  and  source  strength  on  channel  capacity 
is  established  for  the  seismic  signals  and  noise  determined  at  each  site.  We 
find  that  seismic  data  links  of  2  bits/second  are  quite  realizable  for 
distances  of  a  few  kilometers  using  sources  of  less  than  100  Newtons  rms  in 
the  band  2  ^  f  20  Hz. 

The  major  uncertainty  encountered  when  estimating  capacity  is  caused  by 
the  uncertainty  in  Q  and  Che  presence  of  lateral  Inhomogeneltles  even  in  these 
seemingly  siaiple  valleys. 
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3.  SEISMIC  NOISE  IN  BASIN  VALLEYS 

3. 1  Statistical  Classification 

Seismic  noise  is  sensitive  both  to  local  ground  structure  and  Che  kinds 
and  distribution  of  noise  sources.  For  sites  excited  by  a  large  number  of  in¬ 
dependent  sources,  we  can  anticipate  Gaussian  attributes.  If  the  noise 
sources  are  also  well  distributed  about  an  observer,  the  seismic  noise  field 
has  well  defined  spacial  properties  for  uniformly  layered  structures  (6). 

Figure  1  shows  the  distribution  of  the  ground's  particle  velocity 
during  quiet  Intervals  at  a  site  in  Steptoe  Valley  near  Ely,  Nevada.  The  plot 
is  arranged  to  give  a  straight  line  relation  for  Gaussian  variates.  The 

observed  distribution  is  well  represented  by  a  zero  mean,  Gaussian  process 
2  2  2 

N(o,s  )  ,s  ■=  .02  (mlcrons/sec)  (7).  In  remote  areas,  seismic  noise  tends 
Co  exhibit  Gaussian  attributes  (8) .  The  noise  in  the  Basin  and  Range  Province 
sites  is  not  exceptional.  Seismic  noise  in  these  rural  valley  areas  also 
tests  stationary  (7,9).  In  contrast,  noise  in  urban  areas  more  often  is 
non-statlonary,  reflecting  the  rhythms  of  the  work  cycle  (8). 

3.2  Spectral  Estimates 

Being  both  stationary  and  Gaussian,  seismic  noise  at  a  point  can  be 
fully  characterized  by  spectra  (10).  Figure  2  is  a  large  average,  spectral 
estimate  of  base  level  seismic  noise  measured  in  Railroad  Valley.  For 
frequencies  much  above  I  Hz,  the  spectra  are  relatively  flat  when  cast  into 
terms  of  particle  velocity.  The  modest  bulge  in  the  spectra  just  above  10  Hz 
coincides  with  a  loss  of  Impedance  of  the  ground  to  acoustic  loads.  The 
location  of  the  frequency  for  generating  air-coupled  selsmlcs  is  determined  by 
the  velocity  and  density  structure  neighboring  the  site. 
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For  the  sites  considered.  Railroad  Valley  is  much  the  quietest.  The 
seismic  level  here  compares  favorably  to  the  level  forecast  for  quality 
seismic  detection  sites  in  the  southwestern  United  States  (11).  In  turn,  the 
seismic  noise  spectra  at  depth  at  these  sites  can  be  expected  to  be  red 
versions  of  those  measured  at  the  surface  in  that  the  shorter  wavelength, 
higher  frequency  components  are  alAost  always  proportionately  more  attenuated 
at  depth  (12,13). 

A.  SEISMIC  SIGNALS 

A. 1.  Elastic  Response 

For  an  elastic  body  with  homogeneous  boundaries,  the  displacement  at 
time,  T' ,  for  a  point  located  at  X*(x',y',z*)  relates  to  forces  within  a 
source  volume  located  at  X  (x,y,z)  acting  at  time,  T,  through  the  temporal 

convolution  of  the  force  with  the  appropriate  Green  function  (14,15). 

«• 

Uj^{X'  ,T')» ,T';X,T)*f^(X,T)dv  dt 
-  «B  V 

In  elastodynamlcs,  the  convention  is  to  let  a  displacement  at  a  field 
point  caused  by  forces  applied  within  the  source  region  define  response. 

Figure  3.  The  elastic  response  is  causal,  linear  and  time  invariant.  Also, 
when  the  source  region  is  slowly  moving  and  small  with  respect  to  the  radiated 
seismic  velocities  and  wavelengths,  the  source  time  history  can  be  recovered 
by  linear  deconvolution  (16).  Further,  source  and  field  points  can  be 
exchanged  in  accordance  with  the  reciprocity  rule. 

G,^i(X'  ,T';X,T)=  G^j^(X,-T;X',-T’) 

For  uniformly  layered  areas  and  vertical  surface  loads,  ground  response 
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is  Insensitive  to  azimuth.  Placing  the  source  at  the  origin,  we  then  express 
response  solely  as  a  function  of  distance,  r  and  retardation 
time,  t. 


G(r,t) 

4.2.  Normal  Mode  Representation 

We  represent  the  Fourier  transform  of  the  vertical  displace¬ 

ment  of  the  ground  surface  at  a  "large"  distance,  r,  due  to  a  vertical  Impulse 
at  the  origin,  to  be  the  sum  of  the  normal  mode  contributions  attenuated  by  a 

material  loss  factor,  Q, 

N 

G33  (r,u»)=  (oi) -r^exp  (-k^  (a»)  •  r) /2Q  (w)  *  exp  i  (k^  («>)  •  r+j3  ^  («,)  ) 

j-l 

Here  j  is  the  mode  number  and  A,  ^  are  the  apparent  amplitude  and  phase 
associated  with  the  source.  As  given  here,  attenuation  is  governed  by 
cylindrical  spreading  and  a  spacial  Q  loss  defined  by  the  fractional  loss  of 
amplitude  over  one  wavelength.  The  spacial  Q  term  used  here  relates  to  a 
temporal  determination  of  Q  through  the  relation, 

Q~^ (temporal) =  U  .  Q~^ (spacial) 
c 

with  c,  U  the  phase  and  group  velocities  of  the  mode  in  question. 

For  our  test  sites  experimental  values  of  C^^(r,u()  are  fitted  to  a 
fundamental  mode  representation  by  determining  A,k,0,Q  over  a  discrete 
distance  set,  r^,  1,2,3... 
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A. 2.1.  Propagation 

The  average  value  for  unit  amplitude,  k*  weighted  surface  impacts  over 
the  distance,  *  1,2,3,..,L  is 

F33(k'  .«,)=l/L^G33(r^  ,«,)*exp  (-ikr^  )*|  G33(r^,«.)  | 

4^=1 

F33(k',u>)  has  a  maximum  value  of  unity  for  unlmodal,  low  noise  measurements 
when  k'  »  (ui/c).  Figure  4  is  a  sample  computation  of  F33(k',iw)  based  on 
50  response  wavelets  measured  over  the  range  120  ^  r  ^  200  meters.  The 
calculated,  maxima  are  quite  close  to  unity  over  the  pass  band  3.7^f^37  Hz; 
the  seismic  response  at  the  site  can  be  well  represented  by  a  unlmodal  surface 
wave  for  these  frequencies.  Figure  5  locates  the  (k',u))  pairs  that  generate 
absolute  maxima  for  the  data;  the  relation  establishes  the  propagation 
characteristics  at  the  test  site.  The  phase  and  group  delay  times  are 
determined  from  c  °‘ui/k  and  U  ■  dui  /dk 

4.2.2.  Source  Phase 

Once  k'  is  determined,  source  phase  is  calculated  directly  from  the 
residual  phase  at  r>0. 

ARG  ^033(0,  u,  )  ^  =ARG  ^  w  )  ^  •  exp(ik'r) 

In  Figure  6  the  average  phase  residuals  over  reliable  (k,  u>)  pair 
values  are  independent  of  frequency.  The  theoretical  value  for  cyllndrically 
spreading  surface  waves  due  to  a  surface  Impact  on  a  uniformly  layered  elastic 
half  space  at  large  distances,  r,  is  x/4.  The  calculated  source  term  Is 
consistent  with  the  normal  mode  representation  used  here. 
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4.2.3.  Attenuation 


For  unlmodal  propagation  In  uniformly  layered,  lossy  media,  surface 
wave  attenuation  at  a  distance  can  be  obtained  from 

|G(r,tt,)  I  =  A  (ui)  •  1/r^  •  exp(-kr/2Q) 

In  Figure  7  we  plot  In  |  G(r,  50  )r^  |  for  each  of  63  Impact  wavelets 

measured  In  Jackass  Flats  as  they  relate  to  distance.  Q  Is  then  calculated 
from  the  slope  of  the  best  fitting  (least  square)  slope  of  the  data.  A 
measure  of  the  goodness  of  fit,  E  ,  Is  also  calculated  from  the  residuals. 

Figure  8  shows  our  estimate  for  Q  and  E  for  Jackass  Flats  as  they 
relate  to  frequency.  Q  estimates  tend  to  be  decreasing  functions  of 
frequency.  A  Q  of  35  In  surflclal  sediments  Is  modestly  larger  than  Q  values 
used  elsewhere  (1). 

4.2.4.  Prediction  Errors 

The  dispersion  relation  and  Q  values  found  for  Railroad  Valley  are  used 
to  adjust  Individual  Impulse  responses  over  100^ r ^200  meters  to  a  common 
reference  distance,  r>150.  The  average  amplitude  response  for  150  meters, 
[G(150,f)]  ave  Is  then  obtained,  see  Figure  9.  The  ground  Is  found  to  be  a 
band  pass  transmission  element.  The  low  limb  Is  controlled  by  the  nature  of 
the  elastic  response.  The  high  limb  roll  off  Is  dominated  by  the  material 
quality  factor,  Q. 

The  ground  disturbance  predicted  In  Railroad  Valley  for  a  standard 
Impulse  as  It  relates  to  distance  Is  given  In  Figure  10.  The  synthetic 
responses  shown  here  obey  the  A,k,^,Q  parameters  calculated  at  the  site  for 
unlmodal  surface  waves  In  the  range,  I00^r^200  meters. 

The  synthetic  waveforms  are  compared  to  the  actual  wavelets.  Error 


11 


wavelets,  the  difference  between  predicted  and  measured  values,  are  quite 
small.  Figure  11.  The  error  Is  ascribed  to  additive  seismic  noise,  higher 
mode  responses  and  lateral  inhomogenities. 

The  scatter  in  the  individual  amplitude  and  phase  responses  about  an 
average  value  for  Jackass  Flats  is  given  in  Figure  12.  Scatter  diagrams  such 
as  these  are  used  to  support  the  proposition  that  the  original  measurements 
are  heavily  corrupted  by  additive  noise  for  frequencies  greater  than  35  Hz. 
(16). 

It  is  worth  noting  that  the  predicted  selsmics  at  distances  greater 
than  a  few  hundred  meters  are  quite  sensitive  to  our  estimate  of  Q.  A  unit 
difference  in  Q  causes  about  a  6  db  difference  in  the  predicted  midband 
amplitude  response  at  1  km.  Indeed,  a  major  impetus  for  this  study  has  been 
the  lack  of  reliable  attenuation  data  in  Basin  and  Range  Valley  materials. 
Recent  interest  In  mantle  Q  permits  an  upgrading  of  earlier  work  concerning 
cotcmunlcatlng  seismically  at  long  range  (4);  it  bears  little  on  data  rate 
estimates  at  modest  distances  (10km)  for  the  deep  alluvial  valleys  in  the 
Basin  and  Range  Province. 

Even  in  these  supposedly  simple  valley  sites,  attenuation  is  due  to 
lateral  structural  changes  as  well  as  Q  losses.  To  show  this  we  determine  the 
phase  velocity  for  a  number  of  headings  in  Railroad  Valley,  Figure  13.  The 
seismic  response  at  this  site  is  found  to  be  azimuth  sensitive.  Ground 
structure  is  not  laterally  uniform.  Attenuation  is  not  soley  a  function  of 
distance  and  Q  loss.  As  can  be  seen,  the  effect  of  heading  is  somewhat  less 
severe  at  lower  frequencies,  suggesting  the  valley  is  more  uniform  at  depth. 
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5.  CHANNEL  INFORMATION  RATES 

The  capacity  of  a  comntunlcation  channel  In  bics/sec  over  the  band 
a^f^b  is  computed  after  Hartley-Shannon  as 


where  and  are  the  power  spectra  of  the  signal  and  noise  processes, 
respectively  (18).  Channel  capacity  is  a  maxlnnim  when  P^^  «  constant. 

Constant  reception  spectra  can  be  obtained  either  by  weighting  at  the  source 
or  receiver.  In  either  case  we  can  negate  distortions  due  to  the  seismic  path 
by  application  of  a  linear  convolution  operator  that  is  the  inverse  of  the 
ground  response,  G(r,t). 

5.  1 .  Narrowband  Sources 

The  center  frequency  for  minimum  strength  sources  and  optimum 

transmission  is  Illustrated  for  the  transmission  and  noise  characteristics  in 

Railroad  Valley.  Figure  14  Is  the  P  (f)/P  (f)  ratio  for  a  white  source  of 

ss  nn 

strength,  S  =  6.25  Newtons  rms  measured  at  a  range  of  130  meters.  P  /P  is 
®  ss  nn 

a  maximum  near  18  Hz.  In  this  case  a  channel  centered  at  18  Hz  is  optimum  in 
the  sense  that  it  maximizes  channel  capacity  for  a  given  source  level,  band¬ 
width  and  range. 

In  order  that  the  maximum  channel  capacity  for  broadband  sources  be 
realized,  the  Inverse  of  the  path  response  must  be  known.  Even  the  detection 
of  signals  with  a  peak  substantially  less  than  unity  is  not  without 

difficulty.  The  capacity  of  channels,  using  one  or  more  transmitter-receiver 
pairs  uncompensated  for  seismic  dispersion,  passes  a  sub  optimally  coded 
sequence. 
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The  maximum  Information  capacity  for  a  source  of  31.25  Newtons  rms  and 
1  Hz  bandwidth  Is  determined  for  Railroad  valley.  From  Figure  15  we  find  we 
can  transmit  and  receive  In  excess  of  2  blts/secood  over  a  1  Hz  band  at  a 

distance  of  1  km  using  a  source  of  only  31.25  Newtons  rms. 

After  a  range  of  a  few  hundred  meters,  narrowband  sources  are  relative¬ 
ly  efficient  seismic  transmitters.  In  Figure  16  the  effect  of  source  band¬ 
width  and  range  on  Information  capacity  Is  shown  for  a  source  strength  of 
9.375  Newtons  rms.  For  this  source  level  and  locality,  a  bandwidth  of  3  Hz  Is 
optimum  at  a  range  of  0.50  km. 

In  Figure  17  we  determine  the  source  strength  needed  to  maintain  a 

channel  capacity  of  2  bits/second  over  a  bandwidth  of  1  Hz.  A  source  level  of 

at  least  25  Newtons  rms  Is  needed  to  Insure  a  channel  capacity  of  2 
bits/second  at  a  range  of  1  km.  Under  these  constraints  an  approximately 
linear  relationship  exists  between  source  strength  and  distance. 

5.2  Broadband  Sources 

The  Information  rate  for  a  channel  consisting  of  a  broadband  signal 
source  and  a  receiver  whose  response  is  the  inverse  of  that  of  the  ground  Is 
depicted  In  Figure  18  for  a  range  of  150  meters  and  source  strength  of  4.419 
Newtons  In  the  band  of  50  Hz.  The  plot  shows  seismic  communication  to  be  a 
low  pass  phenomenon  In  that  9SZ  of  the  channel  capacity  is  obtained  from  fre¬ 
quencies  less  than  35  Hz.  The  bandwidth  for  efficient  seismic  communication 
Is  also  range  dependent.  To  show  this  we  compute  the  cutoff  frequency  needed 
to  attain  98%  of  channel  capacity  for  white  sources,  Figure  19.  The  useful 
information  bandwidth  for  weak  sources  at  distances  greater  than  I  km  Is  less 
chan  20  Hz  In  Railroad  Valley. 

To  further  show  the  impact  of  range  on  broadband  communications,  we 
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estimate  channel  capacities  In  Railroad  Valley  as  they  relate  to  bandwidth 
over  a  suite  of  distances.  Figure  20.  Broadband  seismic  communications 
require  relatively  strong  sources  to  excite  frequencies  much  above  25  Hz  for 
distances  greater  than  200  meteis. 

6.  TRANSMITTER  AND  RECEIVER  AKRAYS 

Coherent  signals  embedded  In  additive  Incoherent  noise  provide  a 

conservative  basis  to  estimate  the  effect  of  both  transmitter  and  receiver 

arrays  on  channel  capacity.  For  such  a  construction  P  /P  increases 

ss  nn 

directly  as  the  sum  of  the  number  of  receiver  and  transmitter  elements. 

Channel  capacity  for  systems  employing  multiple  transmitter  and  receiver 
elements  Is  shown  in  Figure  21  for  narrowband  transmissions  at  a  range  of  300 
meters. 

The  prospects  of  surface  arrays  separating  signals  from  noise  are 
considerably  better  for  burled  sources  than  for  surface  sources.  The  reason 
for  this  Is  that  both  the  noise  and  signals  excited  by  surface  sources 
propagate  largely  as  low  order  surface  waves.  Arrays  cannot  separate  seismic 
waves  traveling  in  the  same  direction  with  the  same  mode.  In  contrast,  low 
mode  surface  selsmlcs  are  only  poorly  excited  by  buried  sources.  Arrays  using 
velocity  filtering  should  be  quite  effective  for  separating  seismic  signals 
from  noise  in  the  case  of  burled  sources  (19). 

7.  SITE  EFFECTS 

In  Figure  22  the  seismic  wavelets  produced  by  a  standard  Impact  at  a 
distance  of  150  meters  are  shown  for  Jackass  Flats  and  Railroad  valley. 

Seismic  response  Is  clearly  site  sensitive.  It  immediately  follows  that 
channel  data  rates  are  also  site  sensitive.  Optimum  capacity  channels  for  a 
broadband  transmitter  in  Jackass  Flats  is  given  In  Figure  23.  For  the  ranges 
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and  frequencies  considered.  Railroad  Valley  has  the  greater  capacity.  Low 
noise  in  Railroad  Valley  Is  the  controlling  factor. 

All  the  estimates  given  In  this  study  are  for  point  surface  sources  and 
surface  receivers.  For  a  surface  source  and  a  buried  receiver  at  large  offset 
distances  both  the  noise  and  signal  are  attenuated  In  roughly  similar  fashion. 
P  /P  ratios  at  depth  should  be  much  the  same  as  those  found  for  surface 
sources  and  receivers.  The  results  of  this  study  should  apply.  In  contrast, 
Pss/^nn  significantly  smaller  for  a  source  at  depth  and  a  surface 

receiver.  Reliable  assessments  for  communicating  to  the  surface  from  a  buried 
facility  requires  seismic  measurements  at  depth  in  the  areas  of  Interest. 
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FIGURES 


1.  Seismic  Noise  Distribution  -  Steptoe  Valley 

2.  Seismic  Noise  Spectra  -  Railroad  Valley 

3.  Green's  Function  Representation 

4.  Quality  of  Unlmodal  Model 

5.  Measured  Dispersion 

6 .  Source  Phase 

7.  Attenuation:  ETB 

8.  Q  Estimates 

9.  Average  Amplitude  Response 

10.  Predicted  Response  -  Railroad  Valley  (270*) 

1 1 .  Error  Wavelet 

12.  Amplitude  and  Phase  Scatter  -  Jackass  Flats 

13.  Valley  Velocity 

14.  P  /P  Ratio:  RRV  150  Meters 

ss  nn 

15.  Channel  Capacity:  Narrow  Band  Source  RRV 

16.  Effect  of  Bandwidth 

17.  Effect  of  Source  Strength 

18.  Ground  Bandpass  Characteristic:  RRV 

19.  Ground  Cutoff  Frequency 

20.  Channel  Capacities:  Wide  Band  Source  RRV 

21.  Effect  of  Arrays 

22.  Effect  of  Site 

23.  Channel  Capacities:  Wide  Band  Source 
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GLOSSARY  OF  TERMS 


A  Source  Amplitude 

c  Phase  Velocities 

C  Capacity 

db  Decibel 

E  Goodness  of  Fit 

f  Frequency  Hz 

tg  Group  Delay 

tp  Phase  Delay 

G  Green's  Function 

G^^  Vertical  Displacement  of  Ground  Surface  at  a  "Large"  Distance 
1 

j  Mode  Number 

K  Wave  Number 

k*  Trial  Wave  Number 

km  Kilometer 

2 

N(o,s  )  Zero  Mean,  Normal  Process 

P  Power  Spectra  of  Noise  Process 

nn 

P  Power  Spectra  of  Signal  Process 

88 

i 

Q  Material  Quality  Factor 

i  r  Distance 

* 

rms  Root  Mean  Square 

2 

S  Variance 

t  Retardation  Time 

I 

i 
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GLOSSARY  OF  TERMS 


T  Source  Tine 

T’  Time 

t  Wave  Phase  Delay 

P 

tg  Wave  Group  Delay 

U  Group  Velocity 

tu  Angular  Frequency 
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